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Abstract. We present a model for high-energy emission in microquasars where the energy content of the jets is dominated 
by relativistic protons. We also include a primary leptonic component. Particles are accelerated up to relativistic energies in a 
compact region located near the base of the jet, where most of the emission is produced. 

We calculate the production spectrum due to proton and electron synchrotron radiation and photohadronic interactions. The 
target field for proton-photon collisions is provided by the synchrotron radiation in the acceleration region. In models with a 
significant leptonic component, strong internal photon-photon absorption can attenuate the emission spectrum at high energies. 
Depending on the values of the parameters, our model predicts luminosities in the range 10 34 — 10 37 erg s _1 up to GeV 
energies, with a high-energy tail that can extend up to 10 16 eV. In some cases, however, absorption effects can completely 
suppress the emission above 10 GeV, giving rise to different spectral shapes. These results can be tested in the near future by 
observations with instruments like GLAST-Fermi, HESS II and MAGIC II. 
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INTRODUCTION 

Microquasars are binary systems formed by a donor star 
that feeds a compact object (black hole or neutron star), 
where part of the accreted matter is ejected from the sys- 
tem as two collimated non-thermal jets. Depending on 
the mass of the donor star, microquasars are classified 
into high-mass or low-mass systems. The emission spec- 
trum of microquasars covers almost the entire electro- 
magnetic spectrum, from radio wavelengths to hard X- 
rays. Three high-mass X-ray binaries have also been de- 
tected at Tev gamma rays gamma-rays. The gamma rays 
can be produced by interaction of relativistic particles 
in the jet with the radiation field and the matter in the 
winds of the companion star. Several models to explain 
the gamma-ray emission in high-mass microquasars can 
be found in the literature, see for example Ql, |01 and 
|Ht]. However, in the case of low-mass microquasars, the 
density of the radiation and matter fields supplied by the 
donor star is much lower, and the same type of mecha- 
nisms proposed for high-mass systems turns out to be in- 
efficient |4j]. Here we present a model for the high-energy 
emission of low-mass microquasars, based on the inter- 
action of relativistic particles in the jets with the internal 
radiation, matter and magnetic field of the jet itself. The 
model can be also applied to high-mass microquasars, 
taking into account the effect of the external fields. In 
our model the high-energy emission is mostly due to 
hadronic interactions, though we also take into account 
the contribution of primary and secondary electrons. The 



results show that, under certain physical conditions, low- 
mass microquasars can be sources of gamma rays de- 
tectable by the satellite GLAST-Fermi, or Cherenkov 
telescope arrays like HESS II and MAGIC II. 

MODEL DESCRIPTION 

We consider a conical jet, perpendicular to the orbital 
plane of the binary. A scheme of the jet launching region 
is shown in Figure [TJ The jet is launched at a distance 
zo from the compact object and has an initial radius 
ro = O.lzo- The outflow is assumed to be only mildly rel- 
ativistic, with a bulk Lorentz factor Fj et = 1.5. It carries 
a fraction of the accretion power, Lj et = 0.1L a ccr ~ 10 38 
erg s . We further assume that a fraction L m \ = 0.1Lj et 
of the jet kinetic power is in the form of relativistic pro- 
tons and leptons, L rc \ = L p + L e . We relate the injected 
power in protons and electrons through the parameter a, 
Lp — aL e . 

Particles are accelerated by diffusive shock accelera- 
tion leading to a power-law injection, 

Q = Qo— [erg-WV 1 ] . (1) 

z 

The acceleration region is compact, extending from 
the base of the jet up to z max = 5zo- The efficiency 
of the acceleration mechanism is characterized by the 



z * observer 



of both protons and electrons. Both pp and py inter- 
actions create neutral 7T-mesons that then decay to give 
gamma-rays, 
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FIGURE 1. A detail of the jet injection region. The relevant 
geometrical parameters are indicated. 



parameter rj, so that the rate of energy gain for a particle 
of energy E is 



tJi = necB/E. 



(2) 



Here we assume an efficient accelerator with r/ =0.1. 
The magnetic field in the jet decreases as 

B = B -, (3) 

z 

We determine Bo = B(zo) by requiring equipartition be- 
tween magnetic and kinetic energy densities, 
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where vj e t is the jet bulk velocity. Equation yields 
Bo ~ 10 7 G . 

Figure [2] shows the acceleration rate t~l and the cool- 
ing rates = E~ l dE/dt for the different processes of 
energy loss. The only relevant cooling channel for lep- 
tons is synchrotron radiation. This is also true for high- 
energy protons, whereas at low energies proton cooling 
is dominated by inelastic proton-proton (pp) collisions 
and adiabatic losses. The maximum energy of the parti- 
cles is fixed equating acceleration rate and the sum of the 



cooling rates. This condition yields E, 



F SB 10 10 eV 
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10 16 eV and 



Particle distributions in steady state N(E,z) (in units of 
erg - 1 cm~ 3 ) are calculated solving the transport equation 
in the one-zone approximation |5J], 
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(5) 



We consider several photon production mechanisms: 
proton and electron synchrotron radiation, electron rel- 
ativistic Bremsstrahlung, proton-proton (pp) inelastic 
collisions, inverse Compton scattering (IC) and proton- 
photon (py) collisions on the synchrotron radiation fields 
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(6) 



To estimate the spectrum from the n°— decay we fol- 
lowed Refs. [6] and [7] for the cases of pp and p y col- 
lisions, respectively (for more accurate expressions in 
the latter case see [8]). Proton-photon and proton-proton 
collisions also inject high-energy secondary electron- 
positron pairs, product of the decay of charged pions. 
Pairs are also injected directly through photopair produc- 
tion, 
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(7) 



These secondary leptons also contribute to the gamma- 
ray emission through synchrotron radiation. To calcu- 
late the spectrum of pairs, we used the formulae given 
by Refs. |@] and ifToll . The IC spectra were calculated in 
the local approximation of Ref. Ill 111 , whereas for syn- 
chrotron and Bremsstrahlung radiation we used classical 
expressions, see for example Ref. Il2tl . 
All calculations, except those of interactions with mat- 
ter (for these see II 1311 ). were performed in the jet co- 
moving reference frame and the results were transformed 
to the observer frame using the appropriate Doppler fac- 
tor D(0) = [r jet (1 - j3 jet cos ej] "\ We fixed 9 = 30° for 
the viewing angle. See Table[T|for detailed values of the 
model parameters. 

Finally, we studied possible modifications of the pro- 
duction spectrum due to internal absorption effects: high- 
energy radiation in the jet can be self-absorbed by 
photon-photon annihilation, 
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(8) 



We calculated the opacity t(Ej) for a gamma ray to 
escape from the emission region, and then corrected 
the production spectra through the attenuation factor 
exp(— t). The cross section for photon-photon annihila- 
tion and expression for the opacity can be found in 11411 
and 11511 . Cascades are expected to be suppressed by the 
strong magnetic fields, see fltill . 



RESULTS 

Figure [3] shows some of the spectral energy distributions 
(SEDs) obtained for different values of the model param- 
eters. At energies below 1 GeV, the emission is domi- 
nated by the synchrotron radiation of protons and pri- 
mary leptons. The relative importance of the leptonic 
contribution depends sensibly on the parameter a, that 
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FIGURE 2. Acceleration and cooling rates for electrons (left panel) and protons (right panel), calculated at zo for a = 1. 



TABLE 1. Values and ranges considered for the various parameters 
characterizing the jet and the distributions of relativistic particles. 



Parameter 



Value 



Jet injection point 

Jet initial radius 

Size of acceleration region 

Jet bulk Lorentz factor 

Viewing angle 

Jet kinetic power 

Proton-to-lepton energy ratio 

Acceleration efficiency 

Magnetic field at zo 

Minimum proton/electron energy 

Maximum proton/electron energy 



zo = 10 8 cm*_ 

tq — 0. 1 zo cm 

£max ~ 5 ZO 



r, 



jet 



1.5 



= 30° 
Lj et = 1.7 x 10 38 ergs" 
a = 1 - 1000 
7] =0.1 

S = 2 x 10 7 G 



^min(p.e) 
-^maxfp.e) 



2-l00m {pe) c 2 
i 3 x 10 16 /10 10 eVj. 



50Rschw for a 8M., black hole. 
Maximum values attained along the jet. 



fixes the proton-to-lepton total energy ratio. The electron 
synchrotron luminosity ranges from 10 34 — 10 35 erg s _1 
in a proton-dominated jet (cases (a) and (b)), to 10 37 erg 
s _1 in the case of equipartition ( a = 1, case (c)). Case (c) 
is the only case where IC and Bremsstrahlung are signif- 
icant, since the electron population is larger and there is 
an intense synchrotron photon field that serves as target 
for IC scattering. At energies above 1 GeV, in cases (a) 
and (b) the main contribution to the spectrum is due to 
the decay of 7T°-mesons created in pp collisions, reach- 
ing unabsorbed luminosities of up to 10 37 erg s -1 . For 
a = 1, the synchrotron emission of secondary pairs cre- 
ated in py interactions is relevant as well. The target field 
for py collisions is also the synchrotron field of primary 
leptons, and therefore these contributions are not relevant 
for a > 1 . For additional details the reader is referred to 

Absorption effects modify strongly the emission spec- 
trum in the case with a strong synchrotron radiation field. 
In fact, as it can be seen in Figure [3] (d), high-energy 



emission above £y«l GeV is completely suppressed for 
a = 1. In the rest of the cases, the absorption is moderate 
and does not result in significant changes of the emission 
spectrum. 

CONCLUSIONS 

According to the results presented in this work, low- 
mass microquasars could be sources of high-energy ra- 
diation. Currrent atmospheric Cherenkov telescopes are 
not likely to detect them, since they are sensible to pho- 
tons of energy above several hundreds of GeV, and in 
our models emission in this energy range is completely 
absorbed, or the peak in the luminosity is well below the 
detection threshold of the detectors. However, proton mi- 
croquasars could be observed by AGILE and GLAST 
at MeV-GeV energies, and in the future by enhanced 
Cherenkov arrays like MAGIC II and HESS II. Sys- 
tems like those with high-levels of photomeson produc- 
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FIGURE 3. Panels (a) , (b) and (c): production spectra obtained for different values of the proton- to-lepton energy ratio a, 
minimum particle energy £nim> and spectral index of the injection function cc. Panel (d); same spectrum as in (c) corrected for 
absorption. 



tion should also be strong high-energy (>1 TeV) neu- 
trino sources, since neutrinos are not affected by gamma- 
gamma absorption. 
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